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PRINCIPAL  NOTATION 


A Van  Driest  damping  length  parameter 

cf  local  skin-friction  coefficient 

f dimensionless  stream  function 

f u/ue 

K variable  grid  parameter 

L modified  mixing  length,  or  reference  length 

p pressure 

p+  dimensionless  pressure-gradient  parameter 

Rx.  R|_  Reynolds  number,  ugx/v  or  ugL/v 

Rg  Reynolds  number,  uge/v 

u,v  velocity  components  in  the  x-  and  y-directions,  respective! 
uQ  reference  velocity 

u^  friction  velocity  (tw/p) 

x,y  Cartesian  coordinates 

a parameter  in  the  outer  eddy-viscosity  formula,  see  (2.8) 

6 boundary-layer  thickness 

6*  displacement  thickness,  see  (2.5) 

em  eddy  viscosity 

e*  dimensionless  eddy  viscosity 

n transformed  y-coordinate 

transformed  boundary- layer  thickness 

cr> 

e momentum  thickness,  / u/u  (1  - u/u  )dv 

J 0 6 

o 

Ytr  intermittency 

k:  Von  Karman's  constant 


v 


u 

dynamic  viscosity 

V 

kinematic  viscosity 

p 

dens i ty 

T 

shear  stress 

* 

stream  function 

Subscripts 

e 

outer  edge  of  boundary  layer 

w 

wall 

oo 

freestream  conditions 

()• 

differentiation  with  respect  to 

j 


vi 


1.  INTRODUCTION 


Boundary- layer  calculations  are  usually  performed  for  prescribed 
boundary  conditions  which,  for  two-dimensional  incompressible  flows,  may 
be  provided  in  the  form 

u(x,0)  = 0,  v(x,0)  = vw(x),  u(x,6)  = ue(x)  (1.1) 

In  some  problems,  however,  the  external  velocity  distribution  is  unknown  and 
can  be  determined  in  order  to  satisfy  an  alternative  boundary  condition,  such 
as  prescribed  displacement  thickness,  6*(x),  or  the  wall  shear,  tw(x). 

These  solutions  are  obtained  with  "inverse  boundary-layer  methods"  and  find 
application  in  many  practical  problems.  For  example,  the  airfoils  discussed 
in  a recent  article  by  liebeck  (1976)  are  designed  on  the  principle  that, 
in  certain  regions  of  an  airfoil  flow,  the  deceleration  results  in  values  of 
wall  shear  stress  close  to  zero.  Other  important  applications  of  inverse 
boundary-layer  procedures  occur  in  the  calculation  of  duct  flows  such  as  those 
discussed  by  Cebeci  and  Bradshaw  (1977). 

A particularly  important  application  of  inverse  boundary- layer  procedures 
is  their  potential  to  calculate  the  properties  of  separating  and  reattaching 
boundary-layer  flows.  It  is  well  known  that,  for  a prescribed  external 
velocity  distribution,  the  boundary-layer  equations  are  singular  at  separa- 
tion: they  are  not  singular,  however,  when  the  displacement  thickness  is 

prescribed.  This  was  demonstrated  by  Catherall  and  Mangier  (1966)  who  solved 
the  laminar  boundary- layer  equations  in  the  usual  way  until  separation  was 
approached  and  then,  by  assuming  a displacement-thickness  distribution,  cal- 
culated the  external -velocity  distribution  and  local  flow  properties  through 
the  recirculation  region. 

More  recent  investigations  of  the  inverse  boundary-layer  method  have  been 
reported,  for  example,  by  Klineberg  and  Steger  (1974),  Horton  (1974),  Carter 
(1974,1975),  Carter  and  Wornom  (1975)  and  Williams  (1975,1976),  and  have  involved 
solutions  to  the  laminar  boundary-layer  equations  with  prescribed  displacement- 
thickness  and  shear-stress  distributions  in  regions  of  negative  shear  stress. 

The  present  report  is  concerned  with  incompressible  turbulent  boundary- 
layer  flows  with  and  without  separation.  It  represents  a major  extension  to 
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the  work  described  by  Cebeci  (1975a),  which  was  concerned  with  compressible, 
attached  boundary-layer  flows,  and  makes  use  of  Keller’s  (1970)  two-point 
finite-difference  method  to  solve  boundary-layer  equations  appropriate  to 
turbulent  flow.  The  nonlinear  eigenvalue  approach  of  the  earlier  work  was 
found  to  be  inappropriate  to  the  separating  boundary  layers  which  are  thp 
main  concern  of  this  report  and  an  alternative-approach,  based  on  the  solu- 
tion procedures  of  Cebeci  and  Kellei  (1972),  has  been  used.  The  solutions, 
which  relate  to  attached  and  separating  flows,  were  obtained  with  prescribed 
distributions  of  displacement  thickness. 

The  equations  and  boundary  conditions  are  presented  in  the  following  section 
together  with  the  algebraic  eddy-viscosity  formulation.  The  succeeding  sections 
describe  the  solution  method,  referred  to  as  the  Mechul*-function  method,  and 
these  results  and  their  implications.  A brief  statement  of  the  more  important 
conclusions  forms  the  final  section. 


*lit.  unknown  (Turk.) 
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2.  BASIC  EQUATIONS 


2.1  Boundary-Layer  Equations 

The  conservation  equations  appropriate  for  steady  two-dimensional 
incompressible  laminar  and  turbulent  boundary  layers  may  be  written  as 
follows : 


Continuity: 


Momentum: 


3U  i 3 V 
3X  3y 


3u 


u + v ^7 


3U  _ 


3x 


ay 


1 dp  + 1 3t 
p dx  P 3y 


(2.1) 

(2.2) 


t = v |y  - puV  (2.3) 

Corresponding  boundary  conditions,  for  zero  mass  transfer,  are: 

y = ° u,  v = 0 (2.4j 

y h-  cc  u -*■  ug(x) 

With  specified  distributions  p(x)  or  ug(x),  the  equations  and  boundary 
conditions  represent  the  usual  statement  of  the  two-dimensional  boundary-layer 
problem  for  laminar  and  turbulent  flows.  For  convenience,  we  shall  call  it 
the  standard  problem. 


The  inverse  problem,  which  we  consider  here,  is  a flow  with  specified 
displacement  thickness,  <5*(x),  distribution.  Here  6*(x)  is  defined  by 

oo 

6*(x)  = / (l  -~)dy  (2.5) 

0 e 


2.2  Turbulence  Model 

The  solution  of  the  equations  in  section  2.1  requires  a closure  assump- 
tion for  the  Reynolds  shear  stress,  -pu'v'.  In  our  study  we  use  the  eddy- 
viscosity  concept 
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(2.6) 


-pU  V = pf. 


oU 

m ay 


and  prescribe  t in  the  manner  recormended  by  Cebeci  and  Smith  (1974)  and 
shown  by  them  to  represent  various  boundary-layer  conditions  such  as  those  at 
high-  and  low-Reynolds  numbers,  transition  flows,  etc.  According  to  this 
formulation,  the  turbulent  boundary  layer  is  regarded  as  a composite  layer 
consisting  of  inner  and  outer  regions  with  different  eddy-viscosity  formulas 
for  each  region.  In  the  inner  region  the  eddy  viscosity  is  defined  by 

(LmM  L ay  Ytr 


tlere 


A • 26  ^ iT1 


L = 0.4y[l  - exp(-y/A)] 
N = (1  - 11.8p+)1/2  , 


1/2 


tr 


= 1 - exp 


+ vUe  due 
p = 

X_ 

dx 


-G<x  - *tr>  / 

*tr  J 


G = 8.35  x 10"4(u3/v2)  R'1-34 


'tr 


In  the  outer  region  the  eddy  viscosity  is  given  by 


= a 


oo 

/(ue  - u)dy 


' tr 


where 


(2.7) 


(2.8) 


a = 0.0168[1 . 55/(1  + n)] 

II  = 0. 55[1  - exp(-0.243z]/2  - 0.2982^] 


(2.9) 


z]  = R0/425  - 1 
R, 


ue° 
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2.3  Formulation  of  the  Inverse  Problems 


In  our  approach  to  inverse  boundary-layer  problems,  we  use  the  stream 
function  concept  and  write  the  continuity  and  momentum  equations  as  a third- 
order  system  in  dimensionless  quantities  defined  by 


Then  it  can  be  shown  that  (2.1)  and  (2.2)  can  be  written  as 


V: 

v 


(2.10) 


(b*")'  ML-  f (2.11) 

dx  ax  ax 

Here  primes  denote  differentiation  with  respect  to  y and  b is  given  by: 

b = 1 + 4 * 4 = em/v  <2-12) 

The  boundary  conditions  (2.5)  become 


y = 0 


* = s 0 


(2.13a) 


y = ye  *'  - ue  , Ip  = ue[ye  - 6*(x)]  (2.13b) 

where  the  subscript  e denotes  the  boundary-layer  edge  and 

6*(x)  = , u (2.14) 

o 

The  relation  between  ue  and  p is  given  by  Euler's  equation  which,  in  terms 
of  dimensionless  quantities,  can  be  written  as 

du 

^;  + uo  ~ = 0 (2.15) 

dx  e dx 


Two  approaches  to  the  solution  of  the  system  (2.11),  (2.13)  and  (2.15) 
were  attempted.  In  the  first  approach,  which  we  call  the  "nonlinear  eigen- 
value method,"  the  pressure  was  treated,  as  an  eigenvalue  parameter.  The 
second  approach  regarded  the  pressure  as  an  unknown  function;  we  shall  refer 
to  this  method  as  the  "Mechul -function  method."  The  nonlinear  eigenvalue 
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method  is  described  in  detail  in  Cebeci  (1975a)  and  only  a brief  description 
is  presented  here.  The  Mechul-function  method,  which  proved  to  be  more 
appropriate  to  flows  with  separation,  is  described  in  detail  in  section  3. 

For  flows  with  negative  wall  shear,  it  is  necessary  to  make  approxi- 
mations to  the  governing  equations  to  continue  the  calculations  past  the 
separation  point.  Here  we  use  che  approximation  first  suggested  by 
Reyhner  and  Fliigge-Lotz  (1968).  This  approximation  referred  to  as  FLARE 
by  Williams  (1975),  neglects  the  u(3u/3x)  term  in  the  region  of  negative 
u-velocity.  All  inverse  boundary-layer  procedures,  including  the  ones  des- 
cribed here,  use  this  approximation  for  regions  of  separated  flow.  Once  a 
solution  has  been  obtained  with  this  approximation,  however,  the  assumption 
can  be  relaxed  to  allow  solutions  of  the  complete  boundary-layer  equations. 

This  relaxation  of  the  assumption  is  effected  with  a forward  and  backward 
difference  procedure,  as  described  in  section  3.3. 

2.3.1  Nonlinear  Eigenvalue  Method 

Let  us  assume  that  at  x = xQ  we  are  given  the  initial  velocity  profiles 
and  we  seek  a solution  of  the  system  (2.11)  subject  to  (2.13).  In  the  nonlinear 
eigenvalue  method,  Ug  is  assumed  and  with  (2.15)  we  then  obtain  a solution  of 
(2.11)  subject  to  (2.13a)  and  to  the  first  relation  in  (2.13b),  namely  V = ug 
at  y = yg.  Solutions  were  obtained  by  Cebeci  (1975a) for  compressible  flows, 
in  which  the  governing  equations  were  solved  in  transformed  variables  rather 
than  the  physical  variables  considered  here.  For  the  solution  we  can  compute 
6*  (which  we  shall  denote  by  J*)  from  (2.5)  which,  in  terms  of  stream 
function  and  dimensionless  quantities,  can  be  written  as 

«'(*)  ■ <ye  - V“e>  <2J6> 

Recalling  that  the  desired  value  for  the  displacement  thickness  is  s*(x),  we 
form: 

*[ue(x)]  (2.17) 

and  seek  ue  such  that  <t>[ue(x)]  i 0. 


A 
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To  solve  $[up(x)]=  0,  we  use  Newton's  method.  With  an  estimate  u^(x) 
e e 


of  the  external  velocity  we  define  the  sequence  u:  '(x)  by  setting 


— (v+1  ) / v — ( V ) 

Ug  (x)  = Ug  - 


»[Ue(x)3 
34»Lu^(x)  ] /3Ue 


(2.18) 


The  derivative  of  <p  with  respect  to  ue  can  be  obtained  from  (2.17)  by 
making  use  of  the  relation  given  by  (2.16).  This  gives 


ue  \ 


(2.19) 


One  step  of  the  iteration  may  be  summarized  as  follows:  Ug  (x)  is  assumed 

and  a solution  of  the  standard  problem  obtained.  From  this  result  and  from 
the  desired  value  of  ?*(x),  $ is  found  from  (2.17).  The  next  value  of 
ue  can  then  be  calculated  from  (2.18)  provided  that  3<f>/3Ug  is  known.  That 
is,  obtained  by  solving  a set  of  variational  equations  as  described  in  Cebeci 
(1976).  The  iteration  process  is  repeated  until 


^’’(x) -4v)00  <5, 


(2.20) 


where  6-j  is  a small  error  tolerance. 
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3.  MECHUL  FUNCTION  METHOD 


3.1  Numerical  Formulation 


In  the  Mechul  function  method,  we  solve  the  system  (2.11),  (2.13)  with 
the  numerical  method  of  Keller  (1970).  This  is  an  efficient  accurate  two- 
point  finite-difference  method  that  has  been  used  extensively  by  Cebeci  for 
two-  and  three-dimensional  flows  (see,  for  example,  Cebeci  (1974,  1975b)). 
According  to  this  method  we  introduce  new  dependent  variables  u(x,y),  v(x,y) 
so  that  (2.11)  can  be  written  as  a first-order  system.  With  the  function  p 
treated  as  unknown  and  with  the  help  of  the  y-momentum  equation  (-9p/9y  = 0) 
we  can  write  (2.11)  as 


P’  = u 


(bv).  u M _ v li 

V ' 9X  9X  3X 


(3.1a) 

(3.1b) 

(3.1c) 

(3. Id) 


Similarly  (2.13)  becomes 
y = 0 


t = u = 0 


(3.2a) 


y = ye; 


U = ue,  >Pe  = t(e[ye  - /Rj^  <5*(x)] 


(3.2b) 


Here,  for  convenience,  we  have  dropped  the  bars  on  x,  y,  p and  u 


On  the  rectangular  net  shown  in  figure  1,  we  denote  the  node  points  by 
xQ  =0,  x„  = xn  , + kn  n = 1,  2,  ....  N (3.3) 


y0  ■ o, 


’ xn-l  + kn 

yi  = yj-l  * hj 


j = 1*2,  . . . , J j yj  y g 


Nonuniform  net  spacings  kR  and  hj  are  used.  The  quantities  (\p,  u,  v,  p) 
at  points  (xn,  yj)  of  the  net  are  approximated  by  net  functions  denoted  by 
(V?.  u"  v"  p^).  We  also  employ  the  notation  for  points  and  quantities  midway 
between  net  points  and  for  an^  net  function  qj: 


y 


Figure  1.  Net  rectangle  for  difference  approximations. 


xn-1/2  5 \ (*n  + xn-l ^ 

„n-l/2  1 /nn  . n-1 * 

"j  2 7 <1j  * "j  > 


1 

nj-l/2  = 2 


+ 


qj-1/2 


n 

j-1 


) 


(3.4) 


The  finite-difference  forms  of  equations  (3.1a,b,c),  written  for  the  midpoint 
(xn,yj_-|/2)  of  the  segment  P^p2  shown  in  figure  1,  are 


,n  ,n 


"j 

n 

n 

UJ 

Uj-1 

hJ 

n 

n 

Pj 

Jj-l/2 


\H/2 


0 


(3.5a) 


(3.5b) 


(3.5c) 


Similarly  (3. Id)  is  approximated  by  centering  about  the  midpoint  (xn_]/2’ 
y j _ 1 / 2 ) of  the  rectangle  P]P2P3P4  to  obtain 


“n^^v)j-l/2 


v3:1/2^i/2  - (u2)3- 

“ 2an^Pj-l/2  ~ Pj-l/2^ 


1/2 


I;"-1  V0 

P j - 1 / 2 v j - 


- pn_1 

" Rj-1/2 


1/2 


(3.5d) 


o 


Here 


pn-l 

Vl/2  " 


M]-'  - (bvlTj 


'j-1/2- 


(3.6a) 


x„  — x„  , 
n n-1 


(3.6b) 


Similarly  the  boundary  conditions  (3.2)  and  the  relation  (2.15)  become 

♦;  '0.  uj  - u||  - uj,  *"  • i#ye  - «*(x)]  (3.7) 


n <ue>"  n-1 

P + — 2 — = P"  + — - 2 


(3.8) 


3.2  Solution  of  the  Difference  Equations 

If  we  assume  (i>rJ"^ , u1?”^ , v!?"\  p1?"^)  to  be  known  for  0 <_  j ^ J, 
then  (3.5)  for  1 <.  j <,  J-l  and  the  boundary  conditions  (3.7)  yield  a non- 
linear algebraic  system  of  40+4  equations  with  the  same  number  of  unknowns 
(i|ij,  Uj,  vj.  Pj).  To  solve  the  system  we  use  Newton's  method.  We  introduce 
iterates  [ij)(v',  u(v^,  v|v^,  p(v^]  v = 0,  1,  2,  ....  with  initial  values 


(0)  . 


= 0,  = 0, 


(0)  _ 


VW  . yn-l 
o o 


P(S0)  = pjj-1 
0 0 


uj0)  = uj-1,  vj0)  = vj-1,  pj0)  = pj"1  1<j<J-l  (3.9) 

JO)  = n-1  (0)  = n-1  JO)  _ n-1  p(0)  = pn-l  . _ , 

V’j  » UJ  ue  ’ VJ  " VJ  » PJ  \J  J ~ J 


For  the  higher-order  iterates  we  set 

*jV+1)  = ^V)  + ^jV)*  UjV+1)  = ujV)  + a«jV). 

p(v+l)  _ p(v)  + 5p(v) 

j j j 


(v+1)  = (v)  + 6 (v) 


'j 


'j 


j 


(3.10) 


Then  we  insert  these  expressions  in  place 
drop  the  terms  that  are  quadratic  in  [6ijJ 
yields  the  following  linear  system: 


s of  [*,,  u,,  v,,  PJ  in  (3. 
v).  «u  »>,  SvH,  «pSv)].  Th 


(3.5)  and 
This  procedure 


10 


**j  - **j-i  _r  (5Uj  * sum>  * lri>j 

hi 

«uj  - suj-i  -r  (5aj  * * (r3>j-i 


<3.11 ) 


sPj  ~sPj-l  ■ (r4>j-l 


(sl >j«wj  * (s2)jJ,j-l  * (s3>js*j  * ls4>j5*j-l  * (s5>o«uj  * (s6>j«uj-1 

* <s6>j«uj-l  + <s7>J«Pj  * <s7'j«Pj-l  ■ (r2>j 

for  j = 1,  2,  ....  J.  Here  for  convenience  we  have  dropped  the  superscripts 
v and  n,  and  have  reordered  the  equations  so  that  the  AQ  matrix  is  not 
singular.  (For  details,  see  Cebeci  and  Bradshaw,  1977.)  The  coefficients 
(rk)j  (k  = 1 to  4)  are: 

<rl>j  • Vj-1/2  - *J  + *j-l 
<r3>J-l  = hJvj-l/2  " “j  * “J-1 
tr4>j-l  * 0 


(3.12) 


n-1  (bv) j - (bv).  , , 

(r2}j  = Rj-l/2  FTj  + °n  t(<,,v)j-l/2  + Vj-l/2^j-l/2  “ 

“ (lj2)j-l/2  ” ^j-l/2vj-l/2^  “ 2an(Pj-l/2  ~ Pj-l/2) 

The  coefficients  (s^ ) j (k  = 1 to  7)  of  the  differenced  momentum  equation 
Sire: 


(sl}j  = hf  + “?  (,"j  “ ^j'l/2) 

J 


(s2}j  = "bj-l/hj  + T (*j-l  “ ^j-l/2) 


*3)j  = (vj  + vj-l/2] 


(3.13) 


(S4}j  = ~7  (vj-l  + vj-l/2) 
(s5>j  = "°nuj 
^s6^j  ” "anuj-l 


(s7^  = "an 
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The  boundary  conditions  (3.7)  become 


6*0  = h>o  =° 

5u0  - (r2)0  = 0 

(3.14) 

uJ6uJ  + 6PJ  = (r4>J  H 

= (r3)j  -*j  + Uj(yj  -/If  I*) 

(3.14) 

We  use  the  block  elimination  method  discussed  by  Isaacson  and  Keller  (1966) 
to  solve  the  linear  system.  According  to  this  method,  the  linear  system  is 
written  in  the  form 


(3.15) 


(3.16) 
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The  matrices  A . , B j , Cj  are  4x4  matrices  which  are  given  by 


o 

o 

o 

1 

0 0 0 0 

0 10  0 

0 0 0 0 

Ao  s 

o 

C\J 

o 

-C 

1 

1 

o 

III 

o 

0 1 -h0/2  0 

~ 1 

o 

o 

0 0 0 1 

(3.17) 


-1  -hj/2  0 0 

<s4>j  <s6>j  (s2>j  <s7>j 

1 <_j  < J;  Aj  = 

1 -hj/2  0 0 

(s3*j  *s5*j  (sl^j  (s7*j 

0 0 

0 0 

0 -1  -hj/2  0 

0 0 

0 0 

0 0 0 -1 

3.3  FLARE  Approximation  and  DUIT  Procedures 


As  previously  mentioned,  flows  with  negative  wall  shear  require  the  use 
of  the  FLARE  approximation  in  the  region  of  reversed  flow.  According  to  this 
approximation,  and  u^~  are  set  equal  to  zero  in  the  region  of  negative  u-. 
Once  a solution  is  obtained  with  this  approximation,  more  accurate  calculations 
can  be  obtained  without  the  FLARE  approximation  by  using  a forward  and  back- 
ward difference  scheme.  In  the  region  of  reverse  flow  (uj  <^0),  we  write  the 
difference  approximations  of  (3. Id)  for  the  midpoint  (*n+-|/2*  yj-l/2^  the 
rectangle  pi P2P3P4  (see  fl9ure  2)  t0  9et 


(bv)"  - <*)$_, 

FT 


- 3[(*v)j_1/2  + vj-l/2^j-l/2  “ (u<^j-l/2  ^j-l/2Vj-l/2^ 


2*n 


.n+1  .,n 


_2ftrpn+1  - Pn  1 = Rn+I 

j-1/2  j-l/2J  j-1/2 


- Dn+1 


(3.18) 


Here 


Rn+1 

j-1/2 


B * 


(bv)"+1  - (bv)£j  n+1  ? n+1 

FT  + BC(iPv> j_i/2  + (u  ^ j-l/2^ 

J 

1 


xn+l  ~ xn 


■ 


Fiaure  2.  Net  rectangle  for  difference  approximations  for  the  DUIT  procedure. 

With  this  procedure,  the  calculations  are  repeated,  starting  at  the  x-station 
where  the  calculations  first  indicated  reverse  flow  and  are  continued  to  the 
final  x-station,  until  the  solutions  do  not  indicate  any  changes.  This  pro- 
cedure, referred  to  as  DUIT  (downstream  upstream  iteration)  by  Williams  (1975), 
usually  gives  converged  solutions  with  two  or  three  DUITs. 

3.4  Comments  on  the  Solution  Algorithm 

The  solution  of  the  governing  equations  given  by  (2.1)  to  (2.4)  or  (2.11) 
and  (2.13),  requires  initial  conditions  which  can  be  obtained  in  the  trans- 
formed variables. 


n = AL/v'x  y , 4>  = ATvx  f(x,n) 


(3.20) 


Equations  (2.1)  to  (2.4)  become: 


(bf")’  + OLtJ- ff"  + m[l  - (f1)2]  = x f'-~-f"  f£  (3.21) 


T v‘  ' J ~ ' 3X  ■ 3X 

where  primes  indicate  differentiation  with  respect  to  n and 

„ . x due 

* ' car 

e 

Equation  (3.21)  is  subject  to  the  following  boundary  conditions: 


(3.22) 


l 
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n = 0 


f = f'  = 0 


n -*•  K 


f'  - 1 


(3.23) 


Initial  conditions  were  obtained  by  solving  the  governing  equations  for  the 
standard  problem  to  a location  upstream  of  separation. 


The  present  method  has  been  developed  to  allow  nonuni  form  net  spaclngs  in 
the  streamwise  direction  and  across  the  boundary  layer.  In  the  latter  case,  a 
geometric  progression  has  been  used  with  the  property  that  the  ratio  of  lengths 
of  any  two  adjacent  intervals  is  a constant;  that  is,  hj  * Khj_| . The  dis- 
tance to  the  j-th  line  is  given  by  the  following  formula: 

nj  = h1(KJ‘  — 1 )/(K  — 1)  K > 1 (3.24) 


There  are  two  parameters:  h-j , the  length  of  the  first  An  step,  and  K, 
the  ratio  of  two  successive  steps.  The  total  number  of  points  J can  be 
calculated  from  the  expression, 


ln[l  + (K  - lXnyh,)] 

TrTR 


(3.25) 


For  further  details,  see  Cebeci  and  Bradshaw  (1977). 
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4.  RESULTS 


To  test  the  Mechul  function  method  calculations  were  performed  for 
laminar  and  turbulent  boundary  layers  with  and  without  separation,  and 
compared  with  the  results  of  previous  numerical  solutions  for  laminar  flows 
and  with  experimental  data  for  turbulent  flows.  The  latter  comparisons  were 
made  for  a number  of  attached  flows  and  only  for  one  separated  flow. 


4.1  Results  for  Laminar  Flows  with  Separation 

For  laminar  flows  with  separation  we  have  considered  four  separate  flows 
which  allow  the  displacement  thickness  distribution  to  be  specified.  These 
include  the  laminar  separated  flows  of  Williams  (1976),  Carter  (1975)  and 
Briley  (19711.  In  the  last  case,  the  free-stream  boundary  condition  corres- 
oonded  to  a linearly  decreasing  external  velocity  distribution  followed  by  a 
constant  velocity  and  resulted  in  separation  and  reattachment.  Figure  3 allows 
a comparison  of  values  of  local  skin-friction  coefficient  c^  /Rj^  calculated 
with  the  present  method  with  those  obtained  by  Briley  from  the  steady,  two- 
dimensional  form  of  the  Navier-Stokes  equations.  Here  c^  is  defined  by 


c = T*  = 20" (0) 

f 1/2  pu^ 


(4.1) 


and  represents  2ii<"(0).  These  calculations  were  made  for  the  displace- 

ment thickness  distribution  (see  in  figure  3)  deduced  from  the  Navier-Stokes 
solutions.  In  general,  the  agreement  is  very  good  and  the  small  discrepancy 
may  be  associated  with  the  difficulty  of  reading  the  input  7*(x)  distribution 
from  the  graph  presented  by  Briley.  The  Navier-Stokes  solutions  were 
obtained  for  = 2.08  x 10^  and  required  45  minutes  on  a UNIVAC  1108.  The 
results  obtained  by  the  Mechul  function  approach  required  approximately 
10  seconds  on  a CDC-6600  computer. 


We  next  consider  the  two  laminar  flows  with  separation  and  reattachment 
computed  by  Carter  (1975).  These  flows  have  displacement-thickness  distribu- 
tions 6*(x)  given  by 
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0.2 


PRFSENT  METHOD 

+ BRILEY  (1971) 


Figure  3.  Calculated  local  skin-friction  coefficient  distribution  for 
separating  and  reattaching  flow  computed  by  Briley  (1971). 


6*(x)  = 


1 .7208»/x 


2 3 

a1  + a2(x  - x1 ) + a3(x  - x1 ) + a4(x  - x] ) 


o 3 

a,  + a,(x  - Xp ) + a4(x  - x?) 


l.Olxix 


<1  - - , 

Xp  1 x 1 x. 


with 


a-j  = 1 . 7208  ^ 

a2  = (0.5) (1.7208)//^" 

a3  = (0. 5/a-j  ) [6/ A-,  (<5*ax  - a-,)  -4a2] 

a4  = 2/4?  La-,/2  a2  - (<5*ax  - a^]  ^ 

al  " 6max 

V - '/'I^Sax  - *3>  + 

a4  = l/a|  C2(«*,x  - «J)  * V*3'] 

A-j  " x2  x-j  i A2  ~ x3  x2 , X"j  — 1 .065,  x2  — 1 .35,  — 1 .884 

The  first  flow,  referred  to  as  Case  A,  has  = 5.6,  and  the  second  flow 

referred  to  as  Case  B,  has  <5*  = 8.6.  The  two  displacement  thickness 

JTlaX 

distributions  are  shown  in  Figures  5 and  6. 

Comparison  of  present  results  with  those  of  Carter  (1975)  is  shown  in 
Figures  4,  5,  and  6.  The  present  calculations  were  started  at  x = 0 by 
solving  the  governing  equations  in  transformed  variables  for  the  standard 
problem  and  then  at  x = 1 the  Mechul-function  method  was  used  to  solve  the 
inverse  problem  with  the  equations  expressed  in  physical  variables  and  with 
the  FLARE  approximations.  The  procedure  of  improving  the  FLARE  results  by 
the  DUIT  procedure  is  currently  under  examination. 

As  can  be  seen,  the  present  results  agree  well  with  those  of  Carter 
which  made  use  of  the  FLARE  approximations  and  the  DUIT  procedure.  Our 
skin-friction  results  are  in  better  agreement  with  Carter's  results  obtained 
by  the  FLARE  procedure.  His  FLARE  approximation  was  slightly  different  from 
that  considered  here  and  may  account  for  the  disagreement  with  the  present 
results.  The  computer  time  associated  with  the  results  was  again  of  the 
order  of  10s  on  a CDC  6600.  At  almost  all  x-stations  the  convergence  rate 
was  quadratic  and  required  only  two  or  three  iterations,  including  regions 
of  separated  flow. 
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gure  6.  Streamline  pattern  in  separation  bubble  for  Case  B. 


A 


Figure  7 shows  the  results  for  another  laminar  flow,  with  prescribed 
displacement  thickness  distribution,  previously  computed  by  Williams  (1976). 

The  6*-distribution  is  given  by 

6*  = 0.6479(1  + 9 exp  [-0.0625(x  - 14)2]}  (4.4) 

In  this  case  the  flow  starts  at  the  stagnation-point  and  separates  at  approx- 
imately x = 10.  Both  calculations  were  made  by  the  FLARE  approximation.  As 
in  the  previous  cases,  the  solutions  converged  quadratical ly  requiring  only 
two  or  three  iterations,  including  regions  of  separated  flow. 

4.2  Results  for  Turbulent  Flows 

Two  separate  turbulent  flows  have  been  considered.  These  flows  did  not 
have  flow  separation  and  were  computed  with  the  eigenvalue  and  the  Mechul- 
function  methods  to  check  the  predictions  of  both  methods  for  attaching  flows. 
The  two  sets  of  results  were  identical.  The  application  of  the  method  for 
turbulent  flows  with  separation  is  currently  under  examination. 

The  two  attached  turbulent  flows  were  those  measured  by  Schubauer  and 
Spangenberg  (see  Coles  and  Hirst,  1969)  and  labeled  as  44 00  and  4800  at  the 
1968  Stanford  Conference.  Flow  4400  has  a strong  adverse  pressure  gradient 
and  flow  4800  has  a mild  adverse  pressure  gradient.  Figures  8 and  9 present 
the  results  for  these  two  flows.  The  calculations  were  made  first  by  using 
the  standard  procedure:  in  this  case  the  external  velocity  distribution  was 

specified  and  the  boundary-layer  parameters  were  computed.  Next,  the  calcula- 
tions were  made  with  both  of  the  inverse  procedures:  the  measured  displacement 

thickness  distribution  was  specified  and  the  external  velocity  distribution 
computed  along  with  other  boundary-layer  parameters.  As  the  results  show, 
the  agreement  with  experiment  is  very  good;  the  difference  between  calculated 
and  experimental  velocity  distributions  is  small,  and  the  boundary-layer 
parameters  computed  by  standard  and  inverse  procedures  agree  with  each  other. 


1.0  2.0  3.0  4.0  5.0  6.0 


* (ft) 


x (ft) 


Figure  8.  Results  for  flow  4400.  (Inverse  results  were  obtained  by  using 
both  the  nonlinear  eigenvalue  method  and  the  Mechul -function 
method  and  they  are  identical.) 


5.  CONCLUDING  REMARKS 


A numerical  method  for  computing  two-dimensional  Incompressible  laminar 
and  turbulent  boundary  layers  for  separated  flows  is  presented.  The  Reynolds 
shear  stress  term  is  modeled  by  using  the  algebraic  eddy-viscosity  formulas 
developed  by  the  author.  Two  separate  approaches  for  solving  the  governing 
equations  are  investigated.  The  first  approach,  which  we  call  the  "nonlinear 
eigenvalue  method"  treats  the  pressure  as  an  eigenvalue  parameter  and  works 
with  either  physical  or  transformed  variables.  The  second  approach  treats 
the  pressure  as  an  unknown  function  and  seeks  the  solution  of  the  governing 
equations  in  physical  variables.  This  approach,  which  we  call  the  "Mechul 
function  scheme"  works  very  well  for  both  attached  and  separated  flows  while 
the  first  approach  only  works  for  attached  flows. 

Comparisons  and  calculated  results  for  both  laminar  and  turbulent  flows 
show  the  Mechul  function  method  provides  a very  useful  tool  for  computing  the 
separated  flows  with  the  FLARE  approximation.  The  computed  results  agree  well 
with  the  available  numerical  and  experimental  results  and  show  no  signs  of 
numerical  problems  in  regions  of  negative  wall  shear. 
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